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In  this  work,  the  concept  development,  system  layout,  component  simulation  and  the  overall  DOE  system 
optimization  of  a  HT-PEM  fuel  cell  APU  with  a  net  electric  power  output  of  4.5  kW  and  an  onboard  methane 
fuel  processor  are  presented. 

A  highly  integrated  system  layout  has  been  developed  that  enables  fast  startup  within  7.5  min,  a  closed 
system  water  balance  and  high  fuel  processor  efficiencies  of  up  to  85%  due  to  the  recuperation  of  the 
anode  offgas  burner  heat.  The  integration  of  the  system  battery  into  the  load  management  enhances  the 
transient  electric  performance  and  the  maximum  electric  power  output  of  the  APU  system. 

Simulation  models  of  the  carbon  monoxide  influence  on  HT-PEM  cell  voltage,  the  concentration  and 
temperature  profiles  within  the  autothermal  reformer  (ATR)  and  the  CO  conversion  rates  within  the 
watergas  shift  stages  (WGSs)  have  been  developed.  They  enable  the  optimization  of  the  CO  concentration 
in  the  anode  gas  of  the  fuel  cell  in  order  to  achieve  maximum  system  efficiencies  and  an  optimized 
dimensioning  of  the  ATR  and  WGS  reactors. 

Furthermore  a  DOE  optimization  of  the  global  system  parameters  cathode  stoichiometry,  anode  sto¬ 
ichiometry,  air/fuel  ratio  and  steam/carbon  ratio  of  the  fuel  processing  system  has  been  performed  in 
order  to  achieve  maximum  system  efficiencies  for  all  system  operating  points  under  given  boundary 
conditions. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  application  of  compact  fuel  cell  systems  as  range  exten¬ 
ders  for  battery  electric  vehicles  or  auxiliary  power  units  (APUs) 
for  onboard  electric  power  generation  enables  a  significant  reduc¬ 
tion  of  carbon  dioxide  and  other  pollutant  emissions  and  has 
an  improved  NVH  (noise  vibration  and  harshness)  performance 
compared  to  combustion-engine  based  technologies.  Integrated 
fuel  processing  systems  enable  the  utilization  of  the  existing  fuel 
infrastructure  and  have  handling  and  packaging  advantages  over 
hydrogen  based  systems.  The  successful  adoption  of  fuel  cell  APUs 
in  mobile  applications  depends  on  further  technological  develop¬ 
ment  and  cost  optimization  of  the  system  components  and  the 
overall  system.  The  applications  require  fast  startup-times,  a  good 
cyclability,  a  compact  system  design  and  a  minimized  system  com¬ 
plexity  to  reduce  control  effort  and  overall  system  cost. 

Within  this  work,  the  development  of  a  HT-PEM  fuel  cell  APU 
system  with  onboard  fuel  processor  is  presented.  Concept  devel¬ 
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opment,  system  layout,  component  simulation  and  overall  DOE 
(design  of  experiments)  system  optimization  for  maximum  effi¬ 
ciencies  are  described  in  detail.  The  system  is  operated  on  natural 
gas  and  has  an  electric  net  power  output  of  4.5  kW.  Low  pollutant 
and  noise  emissions  as  well  as  a  high  tolerance  against  vibrations 
and  a  high  specific  system  power  density  are  important  criteria  for 
the  application  in  an  inner-city  light  duty  delivery  vehicle. 


2.  HT-PEM  fuel  cell  APU  concept  development 

For  the  APU  development  a  HT-PEM  fuel  cell  with  H3PO4/PBI 
membrane  was  chosen  because  of  its  high  resistance  against  shock 
and  vibration  loads  of  the  mobile  application  as  well  as  shorter 
startup  times  and  improved  thermal  cyclability  compared  to  solid 
oxide  fuel  cells.  Advantages  over  LT-PEM  based  systems  are  the 
increased  CO  tolerance  that  significantly  reduces  the  complexity 
of  the  fuel  processing  system  as  well  as  the  omission  of  anode 
and  cathode  humidification  that,  in  turn,  enables  a  compact  sys¬ 
tem  design.  The  achievable  cell  voltages  of  today’s  state-of-the-art 
HT-PEM  fuel  cells  are  somewhat  lower  compared  to  LT-PEM  and 
solid  oxide  fuel  cells  (compare  Fig.  1 ).  This  can  be  compensated  by 
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Fig.  1.  Current/voltage  characteristics  of  different  LT-PEM,  HT-PEM  and  solid  oxide  fuel  cells  for  hydrogen  operation  and  steam  reforming/autothermal  reformate  operation 
(operating  parameters  and  references  are  listed  in  the  appendix  (Fig.  16)). 


a  reduction  of  the  targeted  current  density  but  results  in  a  tradeoff 
concerning  packaging  volume  and  fuel  cell  cost. 

In  order  to  generate  hydrogen  rich  synthesis  gas  for  the  fuel  cell 
from  a  hydrocarbon  fuel  source,  three  main  technologies  are  used 
for  stationary  and  mobile  applications: 

■  steam  reforming  (SR),  where  hydrocarbons  and  steam  react 
endothermally  to  synthesis  gas 

■  catalytic  partial  oxidation  (CPO),  where  synthesis  gas  is  produced 
via  an  exothermal  partial  combustion  reaction  of  hydrocarbons 
and  oxygen/air  and 

■  autothermal  reforming  ( ATR),  which  is  a  combination  of  exother¬ 
mal  CPO  and  endothermal  SR  and  enables  a  closed  heat  balance  of 
the  reactor  without  any  additional  heat  supplied  to  or  dissipated 
from  the  reactor. 

For  the  APU  application  autothermal  reforming  was  selected 
due  to  packaging  and  transient  performance  advantages  over  steam 
reforming  and  better  durability  as  well  as  reduced  coke  formation 
tendency  compared  to  partial  oxidation. 

On  the  cathode  side  of  the  fuel  cell,  an  increase  in  fuel  cell 
efficiency,  especially  at  high  current  densities,  can  be  achieved  by 
operating  the  system  at  elevated  pressure  levels.  The  higher  power 
densities  that  can  be  realized  for  pressurized  operation  of  the  fuel 
cell  result  in  a  reduction  of  stack  volume,  active  cell  area  and  thus 
fuel  cell  cost.  However,  a  high  compressor  efficiency  and,  possibly, 
also  an  expansion  turbine  are  required  to  minimize  the  required 
compressor  power  input  and  to  enable  a  net  system  efficiency  gain. 
Therefore,  for  APU  systems  in  a  power  range  of  5  kWel ,  low  pressure 
air  supply  systems  operating  close  to  ambient  pressure  are  better 
suited  due  to  lower  system  complexity,  less  control  effort,  reduced 
noise  emissions  and  reduced  air  supply  system  cost. 

In  conclusion,  for  the  development  of  a  fuel  cell  system 
with  an  integrated  fuel  processor  for  mobile  applications  that 
require  intermittent  system  operation,  a  HT-PEM  fuel  cell  with 
an  autothermal  fuel  processor  and  a  low  pressure  air  supply  sys¬ 
tem  enables  the  highest  degree  of  performance  concerning  the 
requirements  system  complexity,  packaging,  degradation,  NVH 
and  startup  time/transient  performance.  It  is  selected  as  the 


basic  concept  for  the  system  development  presented  in  this 
paper. 

2A.  System  layout 

The  system  layout  of  the  HT-PEM  APU  system  is  presented  in 
Fig.  2.  Similar  pressure  losses  of  cathode  and  anode/fuel  processing 
system  facilitate  the  utilization  of  only  one  side  channel  blower  in 
combination  with  an  air  metering  device  and  enable  packaging  as 
well  as  cost  advantages. 

The  autothermal  reformer  is  combined  with  a  2-stage  watergas 
shift  CO  cleanup  unit  enabling  carbon  monoxide  outlet  concentra¬ 
tions  of  <0.7%,  that  result  in  only  a  minor  decrease  of  the  HT-PEM 
performance  at  stack  temperatures  of  160°C.  The  HT-PEM  stack 
is  liquid  cooled  for  fast  and  efficient  startup  and  compact  design. 
It  features  bypass-valves  that  enable  a  disconnection  of  the  anode 
and  cathode  during  startup/shutdown  to  prevent  condensation  of 
water  within  the  fuel  cell.  A  fuel  cell  management  system  has  been 
developed  that  enables  the  measurement  of  cell  voltages  and  pres¬ 
sures/temperatures,  it  monitors  fuel  cell  stack  performance  and 
safety  and  has  outputs  for  system  actuators  like  the  air  supply 
system,  valves  and  relays.  The  combined  startup-/offgas  burner 
preheats  the  fuel  processing  system  and  the  HT-PEM  stack  dur¬ 
ing  startup.  It  converts  the  remaining  hydrogen  in  the  anode  offgas 
during  system  operation.  The  heat  can  be  recuperated  for  the  fuel 
processing  system.  A  low  temperature  coolant  circuit  is  utilized 
to  condense  water  from  the  offgas  burner  exhaust  to  enable  a 
closed  system  water  balance  and  to  cool  the  power  electronics. 
Two  controlled  DC/DC  converters  for  auxiliaries  and  load  enable 
the  integration  of  the  system  battery  into  the  load  management  to 
enhance  transient  performance  and  maximum  power  output  of  the 
APU  system. 

2.2.  Global  system  parameters 

Four  global  system  parameters  can  be  defined  for  the  HT-PEM 
APU  system,  which  affect  system  layout  and  system  operation:  the 
anode  stoichiometry,  the  cathode  stoichiometry,  the  steam/carbon 
ratio  of  the  fuel  processor  and  the  air/fuel  ratio  of  the  fuel  pro¬ 
cessor.  For  all  load  conditions,  those  parameters  can  be  varied 
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LTS:  low  temperature  shift  HTC:  high  temperature  cooling  loop  DC  Load:  DC  converter  load 

Fig.  2.  System  layout  of  the  HT-PEM  fuel  cell  APU. 


within  given  boundaries  (max.  component  temperatures,  required 
heat  flows,  limitations  of  the  chemical  reactions  e.g.  by  the  soot 
limit  or  required  conversion  rates)  to  achieve  maximum  system 
efficiency  (cf.  chapter  4).  The  global  system  parameters  result  in 
massflow,  pressure  loss,  operating  temperature  and  gas  concen¬ 
tration  requirements  for  the  subsystems  and  can  be  defined  as 
follows: 

2.2 A.  Anode  stoichiometry  (pan 

The  anode  stoichiometry  (pan  is  equal  to  the  hydrogen  massflow 
at  the  anode  inlet  mH2,an  divided  by  the  hydrogen  massflow  that 
is  converted  within  the  fuel  cell  stack  rhH2,an,st.-  For  a  given  fuel 
processor  efficiency  i/ref  this  ratio  is  directly  proportional  to  the 
thermal  power  input  into  the  fuel  processor  Pref  divided  by  the  stack 
current  Jsta. 


f^H2,an  f^H2,an  href^ret 

^H2,an,st.  X(/Sta/2F)MH2  ~  /sta 


(1) 


2.2.2.  Cathode  stoichiometry  Xcat 

The  cathode  stoichiometry  Acat  is  equal  to  the  air  massflow  at  the 
cathode  inlet  mair  cat  divided  by  the  air  massflow  that  is  required 
for  hydrogen  conversion  within  the  fuel  cell  stack  mair  cat  st . 

_  mair,cat  _  ^air,cat  ^air,cat 

cat  “  mair,catjSt.  “  x(W4F)(Mo2/?o2,air)  1 


divided  by  the  molar  flow  of  carbon  supplied  to  the  fuel  processor 
via  the  hydrocarbon  fuel  xhCxHy,ref- 


Xref  — 


”H20,ref 

x*icxHy,ref 


^H20,ref 

^CH4,ref 


(4) 


2.3.  Subsystems 
2.3.1.  HT-PEM  fuel  cell 

For  a  given  membrane  electrode  assembly  (MEA)  and  a  defined 
stack  and  flowfield  design,  power  density  as  well  as  cell  voltage, 
thus  stack  efficiency,  are  determined  by  the  operating  parame¬ 
ters:  stack  temperature,  current  density,  CO  concentration  and 
anode-/cathode  stoichiometry.  The  maximum  stack  temperature 
is  defined  by  MEA  degradation  and  thermal  stability  of  the  selected 
stack  materials.  The  current  density  at  the  nominal  operating  point 
depends  on  packaging  and  cost  restrictions  for  the  fuel  cell  stack. 
Carbon  monoxide  concentration  and  anode-/cathode  stoichiom¬ 
etry  have  to  be  optimized  on  system  level  to  maximize  the  net 
APU  efficiency:  cell  efficiency  is  reduced  for  increasing  CO  concen¬ 
trations  due  to  adsorption  of  carbon  monoxide  on  the  platinum 
catalyst.  However,  the  fuel  processor  efficiency  increases  for  higher 
CO-concentrations.  Concentration  losses  of  the  fuel  cell  stack  are 
decreased  at  higher  stoichiometries  of  the  anode  and  cathode 
(Fig.  3),  though  on  system  level  the  power  input  of  the  air  supply 
and  the  fuel  processing  system  increases. 


2.2.3.  Air/fuel  ratio  fuel  processing  system  \ref 

The  air/fuel  ratio  of  the  fuel  processing  system  Xref  is  equal  to  the 
molar  oxygen  flow  at  the  fuel  processor  inlet  ho2,ref  divided  by  the 
molar  oxygen  flow  required  stoichiometrically  for  total  oxidation 
of  the  supplied  hydrocarbon  fuel  h02,ref,st.- 

,  _  ^02,ref  _  ^02,ref  _  ^02,ref 

ref  ^02,ref,st.  (x  +  (y/4))hCxHy,ref  2”CH4,ref 

2.2.4.  Steam/ carbon  ratio  fuel  processing  system  Xref 

The  steam/carbon  ratio  of  the  fuel  processing  system  is  equal 
to  the  molar  flow  of  water  supplied  to  the  fuel  processor  hH2o,ref 


2.3.2.  Fuel  processing  system 

The  fuel  processing  system  consists  of  an  autothermal  reformer 
and  two  watergas  shift  (WGS)  stages  that  reduce  the  carbon 
monoxide  concentration  of  the  synthesis  gas  to  approx.  0.7%.  The 
operating  range  of  the  fuel  processing  system  parameters  Xref  and 
Xref  is  defined  by  the  soot  limit,  the  maximum  catalyst  temperature, 
the  minimum  required  CH4  conversion  rate  and  the  maximum  CO 
outlet  concentration  of  the  ATR  (for  a  given  watergas  shift  layout). 
Operating  boundaries  can  be  determined  with  an  equilibrium  cal¬ 
culation  of  reactor  temperature  and  gas  concentration  by  assuming 
an  adiabatic  reactor  and  a  given  inlet  temperature  (Fig.  4). 
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Fig.  3.  Influence  of  cathode  stoichiometry  on  cell  voltage  (a)  H3PO4/PBI,  1  bar  abs., 
160°C,  reformate/air,  anode  stoich.  1.25,  current  density  0.4  A  cm-2,  CO  0.7%,  own 
data;  (b)  [1];  (c)  [2]. 


Cathode  stoichiometry  /  (-) 

Fig.  5.  Operating  map  of  the  anode  offgas  burner  (complete  oxidation  of  H2  and  CO, 
inlet  temperature  160°C,  Xref=l-5). 


2.3.3.  Startup/offgas  burner 

The  combined  startup/offgas  burner  consists  of  a  catalytic 
burner  for  conversion  of  the  anode  offgas  during  system  opera¬ 
tion  and  a  flame  burner  for  preheating  of  the  system  during  startup. 
With  the  flame  burner,  the  system  can  be  preheated  within  7.5  min. 
During  system  operation,  the  operating  range  of  the  anode  offgas 
burner  is  defined  by  the  minimum  oxygen  concentration  that  is 
required  for  total  oxidation  of  the  anode  offgas  hydrogen  content, 
the  thermal  degradation  limit  of  the  catalyst  at  high  anode  stoi¬ 
chiometries  and  the  minimum  required  temperature  for  methane 
conversion  and  preheating  of  the  reactants.  With  an  increasing 
steam/carbon  ratio  of  the  fuel  processor,  the  operating  range  shifts 
to  higher  anode  ratios  due  to  an  increased  inert  gas  content  in  the 
anode  offgas  (Fig.  5). 


2.3.4.  Air  supply 

The  decrease  in  cell  voltage  due  to  the  power  consumption  of  the 
air  supply  system  can  be  determined  with  the  cathode  stoichiom¬ 
etry  Acat,  the  pressure  ratio  of  the  compressor  77,  the  compressor 
inlet  temperature  T^omp  and  the  overall  compressor  efficiency  rjc omp 
(product  of  isentropic,  mechanical  and  electrical  efficiency). 


AUceu  — 


^o2cpT'omp(nj 


k-1/k 


-1) 


0,  232  x  4Frjcomp 


(5) 


Due  to  similar  pressure  losses  of  cathode  and  anode/fuel  pro¬ 
cessing  system,  one  air  supply  system  is  used.  The  air  metering 


0.00  0.05  0.10  0.15  0.20  0.25  0.30  0.35  0.40 

Air/fuel  ratio  FPS  /  (-) 


device  for  the  anode/fuel  processor  air  massflow  has  to  ensure  a 
high  control  accuracy  of  ±1%  of  the  total  air  massflow  in  order 
to  prevent  fluctuations  of  the  ATR  temperature,  which  strongly 
depends  on  the  air/fuel  ratio  of  the  fuel  processing  system. 

2.3.5.  Cooling  system 

A  low  temperature  coolant  circuit  is  used  in  order  to  cool  the 
power  electronics  and  to  condense  water  from  the  APU  system’s 
exhaust  and  enable  a  closed  system  water  balance.  The  steam  par¬ 
tial  pressure  in  the  anode  offgas  burner  exhaust  depends  on  the 
steam/carbon  ratio  Xref  and  the  air/fuel  ratio  of  the  fuel  processor 
Aref,  the  cathode  stoichiometry  Acat  and  the  stack  current  Jsta: 

2  T  Xref 

P^O  |  +  Xref  +  (2Aref  +  xACat(fsta/4F))(l  +  ( V'’N2,air/V'’02,air))^amb 

(6) 


In  order  to  enable  a  closed  system  water  balance,  the  maximum 
water  content  in  the  APU  exhaust  gas  has  to  match  the  amount  of 
product  water  that  is  generated  within  the  fuel  processor,  the  fuel 
cell  stack  and  the  anode  offgas  burner  during  the  methane  oxida¬ 
tion.  Therefore,  the  maximum  steam  partial  pressure  in  the  APU 
system’s  exhaust  is  independent  from  the  steam/carbon  ratio  of 
the  fuel  processor: 


2 

PH20,exh,max  ^  _|_  (2Xref  +  xAcat(fsta/ 4F))(1  +  (^2, air/ V'02,air))^amb 

(7) 


The  amount  of  heat  that  has  to  be  rejected  from  the  low  tem¬ 
perature  coolant  circuit,  however,  depends  on  the  steam/carbon 
ratio  of  the  fuel  processing  system  and  can  amount  to  up  to  25%  of 
the  thermal  methane  input  power  of  the  APU  system  (Aref  =  0.27; 
fc 0  =  0.5%;  Xref  =  3). 

The  high  temperature  coolant  circuit  is  operated  at  tempera¬ 
tures  of  up  to  1 60  °C  with  heat  transfer  oil  and  is  used  for  preheating 
of  the  stack  during  startup  as  well  as  heat  rejection  from  the  stack 
and  conditioning  of  the  anode  and  cathode  feed  streams  during 
operation. 

2.4.  Influence  of  the  global  system  parameters  on  the  different 
subsystems 


Fig.  4.  Operating  map  of  the  autothermal  fuel  processor  (equilibrium  calculation,  The  influence  of  the  global  system  parameters  on  the  different 

educt  temperature  500  °c).  subsystems  is  listed  in  Table  1 . 
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Table  1 

Influence  of  the  global  system  parameters  anode  stoichiometry,  cathode  stoichiometry,  air/fuel  ratio  fuel  processing  system  and  steam/carbon  ratio  fuel  processing  system 
on  the  subsystems  of  the  APU  (Opt.:  there  exists  an  optimum  within  the  variation  range  of  the  parameter). 


<Pan  t 

^-cat  t 

Xref  t 

^-ref  t 

Input  power  compressor 

t 

t 

t 

t 

Max.  temperature  fuel  processing  system 

** 

t 

Methane  conversion  fuel  processing  system 

Opt. 

t 

Efficiency  fuel  processing  system 

Opt. 

Opt. 

Efficiency  stack 

t 

t 

t 

Temperature/methane  conversion  anode  offgas  burner 

t 

Rejected  heat  HT  coolant  circuit 

t 

** 

t 

Rejected  heat  LT  coolant  circuit 

** 

t 

t 

Fig.  6  shows  the  realized  HT-PEM  fuel  cell  APU  demonstrator 
system  with  integrated  fuel  processor. 


2(H  -  M)k-H4X2M  +  2H+  +  2e“ 

CO  -  M  +  H2Ok^xM  +  C02  +  2H+  +  2e“ 


(10) 

(11) 


3.  Simulation 

Simulation  based  analyses  of  the  HT-PEM  fuel  cell,  the  autother- 
mal  reformer,  the  watergas  shift  stages  and  the  overall  APU  system 
enable  an  optimized  dimensioning  of  individual  components.  Fur¬ 
ther,  the  simulations  help  limit  component  degradation  through  a 
detailed  understanding  of  component  temperature  profiles  and  the 
optimization  of  the  global  system  parameters  to  achieve  maximum 
system  efficiencies  for  given  boundary  conditions. 

3A.  HT-PEM  simulation 

The  influence  of  the  carbon  monoxide  concentration  on  cell  volt¬ 
age  and  efficiency  of  the  HT-PEM  stack  is  a  key  parameter  for  system 
optimization.  Due  to  decreasing  stack  performance  but  increasing 
efficiencies  of  the  fuel  processing  system  for  higher  CO  concentra¬ 
tions,  an  optimum  carbon  monoxide  concentration  exists,  for  which 
maximum  system  efficiencies  are  achieved. 

For  low  temperature  PEM  fuel  cells,  a  number  of  publications 
are  published  that  describe  the  adsorption  of  carbon  monoxide  on 
the  platinum  catalysts  and  the  resulting  decrease  of  free  catalyst 
sites  that  are  available  for  the  hydrogen  oxidation  reaction  [3-7]. 
In  order  to  model  the  influence  of  carbon  monoxide  on  the  HT-PEM 
fuel  cell,  the  approach  of  Baschuk  and  Li  [8]  is  adapted  to  determine 
the  anodic  overpotential  due  to  carbon  monoxide  poisoning  for  the 
HT-PEM  fuel  cell. 

The  reaction  kinetics,  occurring  within  the  anode  catalyst  layer, 
include  adsorption,  desorption  and  electrooxidation  of  hydrogen 
and  carbon  monoxide: 

^H,ads 

H2  +  2M  <=>  2(H  -  M)  (8) 

bH,adsfcH,ads 

^CO,ads 

CO  +  M  ^  CO-M  (9) 

bCO,adsfcCO,ads 


Fuel  processing  system 


Start-/offgasburner 


Side  channel  blower 


Methane  and  Dl-water 
dosing 


voltage 
monitoring 


HT-cooler 


HT-PEM  stack 

System  controller 


For  the  rate  equations,  a  Langmuir  kinetic  is  used  for  the  adsorp¬ 
tion  and  desorption  of  hydrogen: 

rH,ads  =  kH,adscH2(l  ~  ©H  ~  ©cof  ~  ^H,ads^H,ads(^H)2  02) 

For  the  adsorption  and  desorption  of  CO  an  influence  of  the  car¬ 
bon  monoxide  coverage  on  the  adsorption  enthalpy  is  observed  so 
that  the  adsorption  can  be  described  with  a  Temkin  kinetic  [9]: 

ArC0,ads  =  fccO,adsCco(l  ~  ©H  ~  © CO )  exp 

-bco,adskco,ads@co  exp  ^  (13) 

The  kinetic  simulation  parameters  are  adapted  to  test  bench 
results  of  the  HT-PEM  fuel  cell  system  that  is  operated  on  real 
syngas  from  the  fuel  processing  system  (Table  2). 

Reference  measurements  of  the  cell  voltage  for  carbon  monox¬ 
ide  concentrations  of  0.5%  in  the  ATR  syngas  are  compared  with 
voltage  measurements  at  higher  carbon  monoxide  concentrations. 
The  resulting  voltage  difference  between  the  measurements  is 
shown  in  Fig.  7.  The  HT-PEM  simulation  model  enables  a  predic¬ 
tion  of  the  carbon  monoxide  influence  on  the  HT-PEM  voltage  for 
a  wide  range  of  temperatures  and  current  densities. 

3.2.  ATR  simulation 

For  the  concentration  and  temperature  profile  simulation  of 
the  autothermal  reformer  and  the  watergas  shift  stages,  one¬ 
dimensional,  quasi-stationary  reactor  models  are  developed  based 
on  the  conservation  equations  for  mass  and  energy.  The  monolithic 
catalysts  are  subdivided  into  discrete  volume  elements,  for  which 
the  concentration  and  temperature  dependent  reaction  rates  are 
determined. 

For  the  kinetic  simulation  of  the  autothermal  reformer,  the  reac¬ 
tion  mechanisms  for  total  oxidation  (14),  steam  reforming  (15)  and 
watergas  shift  reaction  (16)  are  considered. 


CH4  +  202^C02  +  2H20 

(14) 

CH4  +  H20-e^C0  +  3H2 

(15) 

C0  +  H204^>C02  +  H2 

(16) 

The  approach  ofTrimm  and  Lam  [12]  is  used  for  the  total  oxida¬ 
tion  reaction  of  methane.  For  the  steam  reforming  reaction  and  the 
watergas  shift  reaction  rate  equations  from  Numaguchi  and  Kikuchi 
[13]  are  adapted. 

_  ^toxiPch4Po2  /<tox2Pch4Po2 

rT0X  "  (1  +Kc»Pch4  +^Po2)2  0  +*cVch4  +K8>2) 


^  Py©co  ^ 


Fig.  6.  HT-PEM  APU  demonstrator. 


Table  2 

Kinetic  parameters  of  the  HT-PEM  simulation. 
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Literature  data  [  1 0,1 1  ] 

HT-PEM  simulation,  adapted  to  testbench  data 

H2  adsorption 

/<°H,ads  (ms”1) 

29.8 

29.8 

£H,ads(kJmo1"') 

10.4 

11.7 

H2  desorption 

b°uMs  (xlO11  molirr3) 

4.18 

4.18 

£S,ads  CkJ  mo<_1 ) 

87.9 

87.9 

H2  electrooxidation 

k°H,ox  (mol  ITT2  S”1) 

23.1 

28.1 

£H,ox  (KImo1 

16.7 

10.0 

«h(-) 

0.5 

0.5 

CO  adsorption 

/<°co,ads  (xl04ms-1) 

3.01 

3.01 

£ro,ads(kJmoM) 

47.3 

47.3 

Pi-) 

0.1 

0.1 

y  (kj  mol-1) 

<100°C:  39.7 

135°C:  61.6 

100-115  °C:  41.4 

145  °C:  63.6 

>115°C:  56.5 

150°C:  65.5 

155°C:  67.8 

160°C:  70.6 

CO  desorption 

b°co,ads  (xlO3  molirr3) 

6.87 

6.87 

£CO,ads(kJm01-’) 

100 

100 

CO  electrooxidation 

fc°co  ox  ( x  1 012  mol  m-2  s_1 ) 

3.4 

3.4 

fcO.oxM™1"1) 

127.0 

170.0 

Ole O  (-) 

0.5 

0.5 

rWDRl  =  PCH^DR1  PH2 0  DR1  fcwDRl  (  PCH4  -  7~ - J7 - n  )  (18 ) 

42  \  lPH20^eq,WDRl )  I 

rwcs  =  Ph^o  GS/<  wgs  f  Pc  o  -  („PH2fC°2^  s)  (19) 

2  \  lPH20^eq,WGSj  J 

The  kinetic  parameters  for  the  rate  equations  are  determined 
through  detailed  temperature  measurements  within  the  mono¬ 
lithic  catalyst  (Table  3). 

Fig.  8  shows  the  measured  and  simulated  temperature  profile 
of  the  autothermal  reformer.  At  the  catalyst  inlet,  high  reac¬ 
tion  rates  of  the  exothermal  partial  oxidation  reaction  result  in  a 
steep  increase  in  catalyst  temperature.  The  superposition  of  the 
endothermal  steam  reforming  reaction  with  lower  reaction  rates 


yields  a  temperature  maximum  in  the  foremost  part  of  the  ATR 
catalyst.  With  the  validated  simulation  model,  the  gas  concentra¬ 
tion  profiles  within  the  catalyst  can  be  determined  as  well,  enabling 
an  optimized  catalyst  dimensioning. 

One  of  the  main  influencing  factors  on  catalyst  degradation  is  the 
maximum  catalyst  temperature,  which  depends  on  inlet  temper¬ 
ature,  air/fuel  ratio  and  steam/carbon  ratio  of  the  fuel  processing 
system.  The  simulation  model  enables  an  accurate  prediction  of 
the  maximum  catalyst  temperature  within  ±10I<  for  the  investi¬ 
gated  operating  range  of  the  ATR  (Xref  =  1 .8-3.0,  A,ref  =  0.22-0.30; 
Fig.  9). 

Due  to  the  integration  of  the  offgas  burner  heat  for  educt  pre¬ 
heating  of  the  fuel  processor,  high  conversion  efficiencies  of  up  to 
85%  from  methane  to  hydrogen  are  achieved  with  the  autothermal 
fuel  processor. 


Fig.  7.  Influence  of  CO  concentration  on  the  anodic  overpotential  -  comparison  of  measurement  and  simulation  (H3PO4/PBI  membrane,  1  bar  abs.,  reformate/air). 
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Table  3 

Kinetic  parameters  of  the  ATR  simulation. 


Literature  data  [14,15] 

ATR  simulation,  adapted  to  testbench  data 

Total  oxidation  reaction 

/<°toxi  (xlO11  mol  bar-2  kgcat-1  s-1) 

8.11 

10.11 

£a/roxi  (kjmol-1) 

86.00 

79.50 

k° TOX2  ( X 105  mol  bar-2  kg^r1  s-1 ) 

6.12 

6.82 

Ea,TOX2  (kj  mol-1 ) 

86.00 

79.50 

*8£  (bar-M 

1.26 

1.26 

AH°ch4  (kjmol-1) 

-27.30 

-27.30 

K™’°  (bar-1) 

7.87 

7.87 

AH°o2  (kjmol-1) 

-92.80 

-92.80 

Steam  reforming  reaction 

/<°WDR1  ( x  105  mol  bar-0’404  kg^r1  s-1 ) 

2.58 

2.30 

^a.WDRl  (kjmol-1) 

106.90 

91.00 

«WDR1  (-) 

0 

0 

A/VDR1  (“) 

0.596 

0.596 

Watergas-shift  reaction 

k°wcs  (xlO2  mol  bar-1  kgcat-1  s-1) 

2.45 

84.50 

£a,WGS  (kj  mol-1 ) 

54.30 

30.00 

^WGS  (-) 

0 

0 

Fig.  8.  Temperature  and  concentration  profile  of  the  autothermal  reformer  -  comparison  of  measurement  and  simulation  (ATR  46  mm  x  46  mm  x  1 15  mm,  900  CPSI,  GHSV 
72,000  h-1 ,  A. ref  =  0.3,  xref  =  2.9). 


Fig.  9.  Influence  of  the  air/fuel  ratio  and  the  steam/carbon  ratio  of  the  fuel  processing  system  on  the  maximum  reactor  temperature  -  measurement  and  simulation  (educt 
temperature  450  °C). 
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Fig.  10.  Temperature  influence  on  the  CO  conversion  rate  of  the  low  temperature  shift  -  comparison  measurement  and  simulation  (GHSV  1 0,000  h_1, 1  bar  abs.,  A,ref  =  0.28, 
Xref  =  2.5). 


3.3.  WGS  simulation 

The  reaction  rate  in  the  WGS  reactors,  at  significantly  lower 
temperatures  compared  to  the  autothermal  reformer,  can  be 
approximated  by  the  following  power  law  [16]: 

PcoPH2oP(o2PH2k“  eXP  (j§) 

m  (20) 

Simulation  results  of  the  carbon  monoxide  conversion  within 
the  low  temperature  shift  reactor  are  in  agreement  with  testbench 
data  (Fig.  10).  The  high  temperature  shift  inlet  temperature  profile 
varies  significantly  across  the  catalyst,  so  that  the  CO  conversion 
rates  cannot  be  predicted  within  the  ID  simulation.  Modifications 
of  the  heat  exchanger  upstream  of  the  high  temperature  shift  allow 


for  a  more  homogenous  temperature  distribution  and  a  decrease  of 
the  reactor  volume  in  future  development  stages. 

4.  DOE  optimization  of  the  APU  system 

In  order  to  optimize  the  global  system  parameters  cathode 
stoichiometry,  anode  stoichiometry  as  well  as  air/fuel  ratio  and 
steam/carbon  ratio  of  the  fuel  processing  system  for  all  system 
operating  points,  the  simulation  models  of  ATR,  WGS  and  HT- 
PEM  are  combined  in  an  overall  system  model.  The  influence  of 
additional  operating  parameters  and  the  required  input  power  for 
auxiliary  components  is  incorporated  in  performance  maps  that 
are  based  on  testbench  data.  Design  of  experiments  (DOE)  is  used 
to  generate  cubic  polynomial  models  of  all  system  output  parame¬ 
ters.  The  polynomial  models  enable  the  constrained  optimization  of 
the  system  efficiency  as  a  function  of  net  system  power  for  defined 
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Fig.  11.  DOE  optimization  of  stack  parameters  for  max.  system  efficiency. 
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Fig.  12.  DOE  optimization  of  fuel  processing  system  parameters  for  max.  system  efficiency. 


boundary  conditions  (temperatures,  emissions,  operating  bound¬ 
aries  of  subsystems). 

Figs.  11  and  12  show  the  combinations  of  the  global  sys¬ 
tem  parameters  for  which  maximum  system  efficiencies  are 

achieved  ( . dotted  line).  All  four  global  system  parameters 

are  optimized  simultaneously.  At  high  loads  and,  thus,  high  cur¬ 
rent  densities,  increased  anode  stoichiometries  enable  significant 
efficiency  advantages  of  the  HT-PEM  fuel  cell  stack.  Further¬ 
more,  the  higher  percentage  of  hydrogen  that  is  converted  in 
the  offgas  burner  enables  a  reduction  of  the  required  air/fuel 
ratio  of  the  fuel  processing  system.  In  order  to  realize  low  car¬ 
bon  monoxide  concentrations  at  high  space  velocities  within  the 
watergas  shift  stages,  the  steam/carbon  ratio  is  increased.  The 
cathode  stoichiometry  is  decreased  because  the  moderate  rise 
in  cell  voltage,  with  increasing  cathode  massflow,  is  overcom¬ 
pensated  by  the  higher  power  requirement  of  the  air  supply 
system. 


For  low  system  power  outputs  and  low  current  densities,  the 
influence  of  the  anode  stoichiometry  on  stack  efficiency  decreases, 
and  therefore,  the  highest  system  efficiencies  are  reached  for  low 
anode  stoichiometries.  Due  to  reduced  temperatures  of  the  offgas 
burner,  the  air/fuel  ratio  of  the  fuel  processor  has  to  be  increased 
in  order  to  ensure  adequate  temperatures  and  methane  conversion 
rates  within  the  autothermal  reformer. 

In  addition  to  high  system  efficiencies,  exhaust  emissions  and 
maximum  component  temperatures  can  be  defined  as  constraints 
in  the  optimization  of  the  global  system  parameters.  For  the  FIT- 
PEM  fuel  cell  APU  system,  methane  emissions  increase  at  low 
system  loads  for  maximum  efficiency  operation,  due  to  low  anode 
stoichiometries,  and  thus,  reduced  offgas  burner  temperatures  and 
resulting  low  methane  conversion  rates  (Fig.  13,  dotted  line).  In 
a  second  optimization  step  the  methane  emissions  are  limited 
to  100mgkWhel-1,  resulting  in  the  change  of  the  global  system 
parameters  shown  in  Figs.  1 1  and  12  (  dashed-dotted  line). 
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Fig.  13.  DOE  simulation  methane  emissions  and  max.  ATR  temperature. 
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Fig.  14.  DOE  optimization  of  net  system  efficiency. 


Fig.  15.  Sankey  diagram  of  the  HT-PEM  fuel  cell  system  in  part  load  operation 
(21<WelnetJ. 

Temperature  and  methane  conversion  rate  of  the  offgas  burner 
are  increased  due  to  higher  anode  stoichiometries  and  lower  cath¬ 
ode  stoichiometries.  In  a  third  optimization  step,  the  air/fuel  ratio 
of  the  fuel  processing  system  is  decreased  and  the  steam/carbon 
ratio  is  increased  in  order  to  limit  the  maximum  ATR  tempera¬ 
ture  to  850  °C  and  reduce  catalyst  degradation  (Figs.  11-13, - 

solid  line). 

Fig.  14  shows  the  resulting,  optimized  net  system  efficiency  as 
a  function  of  net  system  electric  power  output  for  all  three  opti¬ 
mization  steps.  Maximum  system  efficiencies  of  25%  are  reached 
with  the  HT-PEM  APU  prototype  for  a  net  system  power  output 
of  2  kW.  For  nominal  load,  a  reduced  fuel  cell  stack  efficiency  at 
high  current  densities  results  in  a  total  system  efficiency  of  20%.  At 
low  partload  system  efficiency  also  decreases  due  to  a  dispropor¬ 
tionate  increase  of  the  auxiliary  components’  power  requirement. 
The  additional  boundary  conditions  regarding  methane  emission- 
/ maximum  ATR  temperature  limitations  have  only  a  minor 
effect  of  less  than  0.5  percentage  points  in  the  lowest  partload 
range. 

A  sankey  diagram  for  part  load  operation  is  shown  in  Fig.  15.  In 
order  to  increase  the  overall  system  efficiency,  the  fuel  cell  stack  has 
the  largest  optimization  potential.  For  the  prototype,  a  cell  voltage 


of  500-550  mV  is  achieved  at  nominal  load,  anode  stoichiometries 
of  >1.25  are  required.  An  increase  in  power  density  and  cell  volt¬ 
age  will  significantly  improve  overall  system  efficiency.  Further 
optimization  potential  can  be  realized  by  increasing  the  fuel  pro¬ 
cessing  system  efficiency  and  reducing  the  power  consumption  of 
the  auxiliaries. 


5.  Conclusion 

In  this  work,  the  concept  development,  system  layout,  com¬ 
ponent  simulation  and  overall  DOE  system  optimization  of  a 
HT-PEM  fuel  cell  APU  with  onboard  fuel  processor  have  been  pre¬ 
sented.  The  HT-PEM  fuel  cell  has  the  benefits  of  an  improved 
thermal  cyclability,  shorter  startup  times  and  a  higher  resistance 
against  vibrations  compared  to  solid  oxide  fuel  cells.  Compared 
to  LT-PEM  fuel  cell  systems  with  onboard  fuel  processor,  the 
system  complexity  can  be  significantly  reduced.  An  ATR  fuel 
processing  system  was  selected  due  to  transient  performance 
advantages  compared  to  steam  reformers  and  reduced  catalyst 
degradation  and  soot  formation  tendency  compared  to  partial  oxi¬ 
dation  reactors.  A  low  pressure  air  supply  system  was  chosen 
because  of  low  noise  emissions,  system  complexity  and  system 
cost. 

A  highly  integrated  system  layout  has  been  developed  that 
enables  fast  startup  within  7.5  min,  a  closed  system  water  bal¬ 
ance  and  high  fuel  processor  efficiencies  of  up  to  85%  (ATR 
methane-to-hydrogen  conversion  efficiency  without  balance  of 
plant  components)  due  to  the  recuperation  of  the  anode  offgas 
burner  heat.  The  integration  of  the  system  battery  into  the  load 
management  enhances  the  transient  electric  performance  and  the 
maximum  electric  power  output  of  the  APU  system. 

Simulation  models  have  been  developed  for  the  HT-PEM  fuel 
cell,  the  autothermal  reformer  and  the  watergas  shift  stages.  They 
enable  the  prediction  of  the  carbon  monoxide  influence  on  cell 
voltage,  a  detailed  analysis  of  the  ATR  reactor  temperature  and 
concentration  profiles  for  optimized  catalyst  dimensioning,  and 
prevention  of  degradation  due  to  catalyst  overheating.  With  the 
WGS  simulation,  the  shift  stage  operating  temperatures  can  be  opti¬ 
mized  in  order  to  achieve  maximum  carbon  monoxide  conversion 
rates. 

Design  of  experiments  has  been  used  to  generate  cubic  poly¬ 
nomial  models  of  all  system  output  parameters.  Based  on  these 
models,  the  four  global  system  parameters  cathode  stoichiometry, 
anode  stoichiometry,  air/fuel  ratio  and  steam/carbon  ratio  of  the 
fuel  processing  system  could  be  optimized  for  all  system  operating 
points.  Boundary  conditions  as  methane  conversion  rates  of  the 
offgas  burner  or  maximum  component  temperatures  were  incor¬ 
porated  in  the  optimization. 

In  order  to  increase  overall  system  efficiency  of  the  HT-PEM  fuel 
cell  APU  prototype,  the  fuel  cell  stack  has  the  largest  optimization 
potential.  An  increase  in  fuel  processing  system  efficiency  and  a 
reduced  power  consumption  of  the  auxiliaries  will  enable  further 
efficiency  gains.  Additionally,  with  a  design  optimization  of  the  sys¬ 
tem  components  a  reduction  of  overall  system  weight  and  shorter 
system  startup  times  due  to  lower  thermal  inertias  of  the  stack  and 
the  fuel  processor  can  be  achieved. 
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Appendix  A. 
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Fig.  16.  Current/voltage  characteristics  of  different  LT-PEM,  HT-PEM  and  solid  oxide  fuel  cells  for  hydrogen  operation  and  ATR/WDR  reformate  operation  (curve  18:  own 
data). 
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Glossary  of  terms 

Symbols 

b:  reaction  rate  constant,  various  (depending  on  reaction  equation) 
c:  concentration,  kmol  m-3 

cp:  specific  heat  capacity  at  constant  pressure,  kj  (kg  I<)-1 
Ea:  activation  energy,  kj 
F:  Faraday’s-constant,  As  mol-1 
i :  current  density,  A  cm-2 
I:  current,  A 

k:  reaction  rate  constant,  various  (depending  on  reaction  equation) 

K:  equilibrium  constant,  - 

m:  massflow,  kgs-1 

M:  molar  mass,  kg  kmol-1 

h:  molar  flow,  kmol  s-1 

p,  p,;  pressure,  partial  pressure,  bar 

Ap;  differential  pressure,  bar 

P:  power,  kW 

r:  reaction  rate,  kmol(s  kgcat)-1 
R:  universal  gas  constant,  kj  (kmol  K)-1 
T:  temperature,  K 
U:  voltage,  V 

x:  number  of  carbon  atoms  within  the  hydrocarbon  molecule,  - 
x:  number  of  cells  in  the  fuel  cell  stack,  - 

y:  number  of  hydrogen  atoms  within  the  hydrocarbon  molecule,  - 
Greek  symbols 
a:  symmetry  factor,  - 

a,  p,  y,  8:  exponents  of  reaction  rate  equation,  - 
ft:  symmetry  factor,  - 

y:  interaction  parameter  CO-adsorption,  kj  mol-1 
A;  difference,  - 
conversion  rate,  - 
p:  efficiency,  - 


pre/:  fuel  processing  system  efficiency,  - 
surface  coverage,  - 
k:  isentropic  exponent,  - 
Xrefi  air/fuel  ratio  fuel  processing  system,  - 
A .cat-  cathode  stoichiometry,  - 
mass  fraction,  kgj  kg-1 
17:  pressure  ratio  air  supply  system,  - 
( pan •'  anode  stoichiometry,  - 

Xref •'  steam/carbon  ratio  fuel  processing  system,  molH2o  mole-1 
i/q:  molar  fraction,  molj  mol-1 

Indices 

a:  activation 
ads:  adsorption 
amb:  ambient 
an:  anode 
comp:  compressor 
exh:  exhaust 
FB:  flame  burner 
i:  component  i 
cat:  catalyst 
cath:  cathode 
m:  molar 
max:  maximum 
min:  minimum 
ox:  oxidation 
Ref:  reformer 
SR:  steam  reforming 
st:  stoichiometric 
Sta:  stack 

TOX:  total  oxidation 
WGS:  watergas-shift 

Abbreviations 

APU:  auxiliary  power  unit 

ATR:  autothermal  reforming 

CB:  catalytic  burner 

CPSI:  cells  per  square  inch 

DOE:  design  of  experiments 

FPS:  fuel  processing  system 

GHSV:  gas  hourly  space  velocity 

HT:  high  temperature 

HTS:  high  temperature  shift 

LT:  low  temperature 

LTS:  low  temperature  shift 

MEA:  membrane  electrode  assembly 

NVH:  noise,  vibration  and  harshness 

0GB:  offgas  burner 

PEM:  proton  exchange  membrane 

CPO:  catalytic  partial  oxidation 

SOFC:  solid  oxide  fuel  cell 

SR:  steam  reforming 


